Mixing of scalar glueballs and 
flavour- singlet scalar mesons 

UKQCD Collaboration 

C. McNeile and C. Michael 
Theoretical Physics Division, Dept. of Mathematical Sciences, 
University of Liverpool, Liverpool L69 3BX, UK 

Abstract 

We discuss in detail the extraction of hadronic mixing strengths from 
lattice studies. We apply this to the mixing of a scalar glueball and a 
scalar meson in the quenched approximation. We also measure correla- 
tions appropriate for flavour-singlet scalar mesons using dynamical quark 
configurations from UKQCD. This enables us to compare the results from 
the quenched study of the mixing with the direct determination of the 
mixed spectrum. Improved methods of evaluating the disconnected quark 
diagrams are also presented. 

1 Introduction 

Lattice techniques are well developed to describe mass spectra. What is much 
less well studied are hadronic transitions. Here we shall concentrate on purely 
hadronic transitions such as glueball mixing with scalar mesons, string breaking, 
flavour singlet pseudoscalar mass generation, etc. In full QCD studies on a 
lattice, much as in experiment, one will obtain the mass values of the resulting 
mixed states. By varying quark masses and the number of quark flavours, one 
may be able to go beyond experiment and so help to substantiate or vitiate 
phenomenological models. 

Within a quenched or partially quenched lattice approach, one can in prin- 
ciple learn much more: obtaining estimates of the mixing strengths themselves. 
This is the approach that we analyse in detail. We then apply it to the mixing 
of a glueball with a scalar meson. This is of considerable phenomenological 
interest: the fate of the glueball is widely debated [Q, ||. 

As a counterpoint to our quenched study of this mixing, we also determine 
the mixed spectrum directly for two flavours of degenerate sea quark. This 
provides a check on our approach and, incidentally, indicates evidence for a 
surprisingly light scalar state at the lattice spacing we employ. 

We include an appendix giving details of the variance reduction technique 
relevant to determining disconnected fermionic loops which are needed in our 
study of flavour singlet mesons. 
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2 Lattice Analysis 



Here we discuss the formalism on a lattice to extract hadronic mixing. To set 
the scene, the variational approach is first summarised and the simpler case of 
weak/electromagnetic matrix elements is reviewed. Then we discuss hadronic 
mixing matrix elements. 

2.1 Variational methods 

Consider a hadronic correlator Cij (t) where t is the lattice separation in the time 
direction and i, j label the type of operator used to create/destroy the hadron 
(eg. whether local, fuzzed, etc). We assume there are N types of operator. 

In an ideal world of infinite statistics the matrix Cij(t) of correlations can 
be written in terms of the eigenstates of the transfer matrix as 

C(t) = A T e- mt A (1) 

where the intermediate state sums are over all "particles" allowed. 

However in practice, we only have a small matrix C and because of noise 
there is a truncation of the sum in the above equation. 

One standard approach to circumvent these problems is the variational 
method. This can be motivated by maximising 

UiCij(t) Uj (2) 

subject to constant Uidj(t — l)uj which leads to the generalised eigenvalue 
equation 

C ij (t)uf = X a C ij (t-l)uJ (3) 

with a = 0, ... TV — 1 and where the ground state with alpha = has the 
largest eigenvalue A . It is usual to relate these eigenvalues to masses (energies 
in general) by X a = e~ m =. Given exact data, N eigenvalues can be determined. 

In order to isolate a particular state, usually the ground state, one can use 
these variational eigenvectors to form a new basis. This is 

C af3 (t) = ufCijit)^ (4) 

Then at times t and t — 1, C a p will be diagonal and the diagonal elements will 
decrease like A" (ie as e~" la ) as time increases from t — 1 to t. 

A typical use of the variational method to extract the ground state signal 
is to form this variational basis using a low value of t where statistical errors 
are relatively small and then to explore the t dependence of Coo(i) at larger 
t- values to extract the asymptotic behaviour. Note that the finite (N state) 
variational basis derived at finite t will not match exactly to the true spectrum 
of excited states and hence Coo will have some small remaining contamination 
of excited states. Indeed it is possible to use the variational estimate of the 
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mass gap to the first excited state (mi — mo) to control this extrapolation to 
large t to determine mo. 

In order to make contact with the alternative procedure of fitting C directly 
to M states over some range of t: 



M-l 



= X! c ? c j e ~ mat (5) 

we note that u corresponds to the right eigenvectors of the non-symmetric matrix 
C _1 (t — l)C{t) and one can introduce left eigenvectors v which satisfy 

< = a,-(t-i)u« (6) 

and are orthogonal to the right eigenvectors 



a 



5 a p (7) 



In terms of these eigenvectors, we have 

N-i 



Cii{t)= 1 £v?vf\ a (8) 

and 



N-l 



1) (9) 

Thus we see that the variational method corresponds to making an exact fit to 
the data at t and t — 1 with N states with the eigenvalues giving the masses 
and the left eigenvectors v are the couplings. 



2.2 Matrix elements: Operator insertions 

Here we consider first the simpler case where a three point function is evaluated 
with an explicit operator (current) at an intermediate time t\. This is the 
case of weak and electromagnetic current insertions and also for some hadronic 
studies. One example is the study of semi-leptonic decays fl. In general the 
quantum numbers of the states propagating before and after the insertion may 
be different. We shall assume that the spectrum of states is discrete. In a 
finite spatial volume, the two particle spectrum will indeed be discrete and it 
is possible to make use of this to explore relevant matrix elements Q for two 
particle systems. 

We study the general behaviour of the three point function with insertion at 
t\ where t is the lattice separation in the time direction and i, j label the type 
of operator used to create/destroy the hadronic state at times and t where 
there are N, N' types of operator respectively. We include M states with 
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masses (energies) m a from time to t\ and M' states with masses (energies) m'p 
from t\ to t: 

M-l M'-l 

X ij (t 1 ,t)= ]T E cfe-^^^e-™^*-* 1 )^ (10) 

Q=0 0=0 

The task is usually to determine the matrix element xoo corresponding to the 
ground state hadrons. One can employ a fit to the above expression along with 
fitting the two point functions (C with couplings c and masses m and C with 
couplings d and masses m'). Care should be taken to include sufficient states 
M and M' . A sensible criterion is that M should be chosen so that a good fit 
is obtained to the two point correlator over time interval t\ for the appropriate 
hadron with creation/destruction operators as used in the three point function 
for that hadron and M' should be chosen likewise so that a good fit was obtained 
for the appropriate two point correlator over time interval t — t\. 

Another way to approach this analysis is to use of the variational technique 
discussed above. Then with a variational basis for the appropriate two point 
functions (of size N, N' respectively) one can form 

X a/3 (t 1 ,t) = ufX ij {t 1 ,t)uf (11) 

This expression will have non-zero off-diagonal elements in general since the 
operator insertion need not be diagonal in the spectrum basis. In this variational 
approach, one can extract xoo by taking a ratio to the variational ground state 
two point functions determined as above 

xoo = X QO (t 1 ,t)/yJc o{t)C' 00 {t) (12) 

where this ratio should be independent of t\ provided t\ and t — t\ are not too 
small. 

2.3 Matrix elements: Hadronic mixing 

The problem of determining hadronic matrix elements involved in mixing - for 
example in mixing of glucball and scalar meson, in string breaking or in flavour- 
singlet meson studies - is much less straightforward. In full QCD, mass eigen- 
states can be determined directly and one is able, much as in experiment, to 
determine the masses of the resulting mixed states. In quenched, or partially 
quenched, studies it is possible to study mixing more directly by evaluating 
correlators between the different states involved. This is the area we explore 
here. 

In lattice QCD with Euclidean time, the main factor is that the lightest state 
allowed will dominate the correlation at large time separation t. Thus in a study 
of glueball decay to two pions, for example, the two pion state will be lightest 
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and dominate. In general this makes hadronic decays difficult to study [pi. As 
we shall see, the way forward is to restrict study to transitions in which the 
energies are similar: on-shell transitions. In the case of glueball decay to two 
pseudoscalar mesons, this would imply varying the quark mass until the two 
pseudoscalar mesons had comparable energy to the glueball, as was arranged in 
a pioneering study of this |(| . 

Here we have in mind several situations where energies are comparable: (i) 
scalar glueball mixing with scalar mesons, (ii) B B mixing with a static QQ at 
separation R, (iii) mixing of ss pseudoscalar mesons with singlet qq pseudoscalar 
mesons. In each case one system is fully treated (i.e. the gluonic interactions in 
the scalar glueball and the static QQ or the qq state with sea quarks q in the 
vacuum) while the other has valence quarks which are treated in a quenched (or 
partially quenched) approach. 

Our notation is that we consider M states with masses m a to be fully de- 
scribed (that is to say they are either fermionic with those quarks present in 
the vacuum or they are fully described in a quenched theory as are glueballs or 
potentials). We also consider M' states with masses sp which are quenched or 
partially quenched (i.e. fermionic with valence quarks which do not participate 
in the sea). Here we assume that the spectrum of states is discrete in each case. 
To be more specific, in one application one can think of the m states as glueballs 
and the s states as scalar mesons. Some relevant correlations are illustrated in 

O-O 0--O 

Figure 1: We illustrate correlations among scalar glueballs (created by a closed 
Wilson loop) and scalar mesons made from quark-antiquark. Clockwise from 
the top left: the disconnected fermionic correlation (D 3S ), the cross correlation 
of a fermion loop with a Wilson loop (C ms or H), the correlation of Wilson 
loops (C mm ) and the connected fermionic correlation (C ss ). 

Thus the s states do not contribute to correlators unless an explicit operator 
creates or destroys them whereas the m states can occur as intermediate states 
in any correlator with the correct quantum numbers. We shall treat the s states 
as those given by the relevant connected fermionic correlator. The disconnected 
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contribution to any correlator will be included explicitly. The model used for 
these disconnected contributions may depend on the application. Assuming 
that the disconnected fermionic loops in the s state to s state correlator are 
joined by gauge links, one would expect transitions between these s states and 
m states to be the relevant description. In some cases, however, such as pseu- 
doscalar mesons, the gluonic links between the disconnected loops are expected 
to be short ranged and so are treated as local. In this case an explicit mixing 
coefficient between disconnected s states should be introduced. 

We will assume that transitions between states of mass m a and sp are local 
and have strength x a p . The goal is usually to determine the transition strength 
between ground states: x o- Because the transition is hadronic, there is no ex- 
plicit matrix element insertion. One must deduce the strength of the transition 
from a study of the two point correlators alone. 

Consider now the two point correlators at separation t between operators 
creating either of these types of state. For each type of creation or annihila- 
tion operator, we consider several different operators labelled i,j = 1,N and 
k, I = 1, N' respectively (here we have in mind different smearing or fuzzing pre- 
scriptions). Then we have a description in terms of transfer matrix eigenstates 
as 

M-l 

c™ m (t) - E c ? e ~ mQ S a ( 13 ) 

t M-lM'-l 

=EEE dte-^Xape-'rt-VdZ (14) 

t 1= a=0 0=0 

There will also be additional terms with t — > T — t for a lattice periodic in time 
with extent T - these we do not write explicitly. Now for the quenched s correla- 
tion one can separate it into connected and disconnected fermion contributions 

CU{t)=CU{t) + D%{t) (15) 

with 

M'-l 

E dfe-^df (16) 

0=0 

and 

t M'-1M'-1 t M-l 

0J?(*) = E E E E E^ e " S/3tl ^ e " roo(t2 " tl) ^' e "^ ( *" t2)d f 

ti=0 0=0 P'=0 t 2 =ti a=0 

(17) 

One might also include contributions to D ss coming from a direct (i.e. not via 
an m state) transition at t\ from state sp to spi with mixing strength ypp/. 
As discussed above, this is appropriate for a discussion of pseudoscalar meson 
mixing, for example. In general, if propagation of the states of mass m a over 
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zero time steps is included, then the above formula does include contributions 
giving the effect of such direct transitions. 

The main problem with extracting the mixing matrix elements xoo is in 
removing the excited state contributions. Unlike in the case of the explicit 
insertion at t\, here we have no dependence on t\ (and £2) in the observables. 
This implies that t\ values from to t will be allowed so that there will be 
excited state contributions to the correlator which are not suppressed (since 
e -(mi-m )ti w ith ti = is large, for example). 

The cleanest way to circumvent this is in the case when the masses are 
equal, i.e. too = sq- In this case, called the on-shell case, the mixing observable 
C" k ls (t) will have a contribution from the sum over t\ of c°xooie" m °*d° from 
the ground states whereas the dominant excited state contributions behave as 
c®xoie~ mot dl and will be suppressed by a relative factor of 1/t. This is a much 
smaller relative suppression than the factor of e -( m i-" l o)* w hich applies to two 
point correlators, but it is sufficient to remove excited state contributions. This 
implies that xqq can in principle be extracted unambiguously. Likewise Df,j(t) 
has a contribution of d°xoo£ 2 e~ mo *xood°/2 from the ground states which also 
dominates, by a factor of i, any excited state contributions. This gives a further 
cross-check since Xqo can be determined in these two independent ways. 

The analysis when too 7^ so is more tricky. The only ^-dependences that can 
be observed will be of the form e~" l °* and e~ s ' jt (and also te~ Sf3t for D ss (t)). 
Then assuming that the M, M' state fits to C™ m (t) and C|f (*) yield the masses 

m a , sp and couplings cf and d k of both ground states and excited states, the t 
dependence of and are available to determine x a p- In principle there 
are enough such independent ^dependences to determine the mixing parame- 
ters, given sufficiently precise data. 

One way to see this is to use the variational formalism with eigenvectors uf 
obtained from C™ m (i) and wf obtained from C k f(t) using time values of t and 
t — 1. Then we can project into this variational basis 

H afj {t) = u?C% s (t)w k (18) 

If the variational basis corresponds to the exact spectrum then we would have, 
using continuum evaluation of the sum over t\ as would be appropriate for a 
very small lattice spacing, 

H afj (t) = x afj (19) 

sp - m a 

from which x a p can be extracted. In general, however, the variational basis does 
not correspond to the exact spectrum. A fit to the t-dependence then is needed. 
Provided enough operators TV and TV' are used (namely N > M, N' > M'), 
there is sufficient information to extract the parameters x Q( g in principle. A 
similar variational analysis of D s k \ (t) is also possible and this gives another way 
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to determine constraints on the parameters x a p. lf/3 = /3' then, using continuum 
evaluation of the sums over t\ and t%, 



(20) 



D„r{t) = \^x M ( S0 -m a r 

while if 0^0', then 

E( e -m a t 
Xa/3Xaf3> 7 w 7 
\(sf3 - m a ){s [3 , - m a ) 

+ (s p - m a )(s0 - sp>) + (sp> - m a ){sp - sp)) ^ 

Note that this approach assumes that an accurate description of the diagonal 
two-point correlators exists which is valid down to t = 0. This is in practice hard 
to achieve. In particular the t = correlator is often quite different (sometimes 
even in sign) from the t > correlators. We now address the possibility of 
excited states that contribute only to t = since these are needed to cope with 
this data. 

Let us explore this situation with m ^ So in a simple example. We will 
assume precise data are available for the correlations at all t values and that, 
by an optimal choice of operators, the diagonal correlations C mm and C ss are 
described exactly by one state for t > 1 and so have an additional excited state 
contribution of the form 

C mm (t) = c°e- mot c + cMtoc 1 (22) 

C ss (t) = d°e- sat d° + d^tod 1 (23) 

Here we suppress the operator labels (i,j,k etc) since we are considering the 
case that an optimum combination of them has already been taken to isolate 
the ground state (i.e. just taking a — and = above). Then the cross 
correlation has the form for t > 0, including the excited state contributions 
from t\ = and t\ = t, 

H(t) = J2 c°e- motl a;ooe- so(t - tl) d + c^we'^d + c° 'e'™ * x m d l (24) 
ti 

Now by completing the sum over t\ as a discrete sum, H(t) can be expressed as 

H = c a d°e- m « t (A + BX t ) (25) 

with A = e -(«o-m ) ) 

,4 = ^ + ^ (26 , 
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and 

s = "T^a + ^ (27) 

Here d°, d 1 , c°, c 1 and A are known in principle but the mixing parameters 
Too, #01 an d %io are to be determined. With perfectly precise data for H, only 
the coefficients A and B can be determined, hence the three mixing parameters 
cannot be independently determined. Thus t o cannot be determined, even 
in principle. The exception to this is when A — > 1, since the contribution of 
.Too is then greatly enhanced. In fact in this case Too can be read off from the 
coefficient of the linear term in t in H as discussed above. 

Now when data are available for the correlation D then additional constraints 
exist. For our example this is 

D{t) = d d°e- mot (X + FA* + Zt\*) (28) 

with 

X = (^) 2 + (d^/dy + 2 f^0 (29) 

Y = -A(2 - A)(^o_) 2 - 2 ( 7° X ° A 1 ) y + 2 XwXn d 1 /<P + x 2 w (30) 

and 

Z=-\x 2 m /(l-\)+x 2 w (31) 

Here x\\ is a new mixing: between the excited states in both sectors. We now 
have three additional constraints (X, Y and Z) given accurate data for D(t), 
with only one additional parameter. Thus in the case of this simple model, the 
measurable quantities overdetermine the mixing parameters. Again as A — > 1, 
this expression simplifies and the coefficient of t 2 in D{t) gives Xq directly. 

One way to check that ground state contributions dominate is to extract 
Too from H and D at several t values neglecting excited states and check for 
consistency. So the relevant expressions will be 

x H (t) = V ' =- 32 

y/C mm (t)C as (t) 1 + A + . . . + A* 

and 



x D (t) = J P® - A</2 = (33) 

W \/c-(t) v /l + 2A+... + (i+l)A t V ' 

In our application below A = 0.64 for strange quarks, so the enhancement of 
the ground state mixing (too) by factors of 1/(1 — A) is not very big. However, 
we use information from both H and D which does provide a cross check in 
principle of our assumption that the excited state contributions (such as xoi 
and tio in the above example) are negligible. 
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L a 


K 


ro/a 






12 


0.1390 


3.05 


0.707(5) 


0.78 


1(5 


0.1390 


3.03 


0.701(6) 


0.78 


12 


0.1395 


3.44 


0.558(8) 


0.71 


16 


0.1395 


3.44 


0.564(4) 


0.72 


12 


0.1398 


3.65 


0.476(16) 


0.67 


16 


0.1398 


3.65 


0.468(5) 


0.67 


12 


0.13843 


2.94 


0.736(2) 


0.78 


12 


0.14077 


2.94 


0.529(2) 


0.65 



Table 1: Properties of the N f = 2 lattices from UKQCD fjj with f3 = 5.2 
and Csw — 1-76 and (last two rows) Nf — lattices || with (3 = 5.7 and 
C SW = 1-57. 

In summary, when we have an on-shell hadronic transition we can extract the 
mixing matrix element with a power suppression in t of excited states compared 
to the exponential suppression which applies to extracting masses. When the 
masses are not equal, excited states can still be removed in principle if a precise 
study is made of both correlations of type H and D. Note that in practice, if 
the difference in masses is significant, the extraction of the mixing strength will 
be very difficult. 

3 Lattices 

We choose to study sea quark effects using the configurations with Nf — 2 
at f3 = 5.2 with C SW = 1-76 from UKQCD §. Two spatial lattice sizes are 
available (12 3 and 16 3 ) so that finite size effects can be explored. We use the 
three lightest sea quark masses available and we use valence quark masses equal 
to the sea-quark mass. The lattice information is summarised in Table 1. Local 
and non-local meson operators were used with fuzzing radius 2 with 5 iterative 
levels with coefficient 2.5. 

We also use for comparison a quenched lattice at (3 = 5.7 of size 12 3 24 with 
valence fermions of k — 0.14077 (approximately strange mass)and k = 0.13843 
(approximately twice strange mass) with Csw = 1-57, as studied previously 

In order to improve the statistics we measure the disconnected diagrams on 
configurations separated by less trajectories than for the connected correlators as 
shown in Table 2. Even though there may be some autocorrelation among these 
measurements separated by less trajectories, we find that this approach does 
allow the statistical error to be reduced. Indeed this is the approach that was 
used in glueball studies, where the measurement time is very small so one might 
as well measure almost every configuration - indeed we follow this approach here 
when considering glueballs, as discussed later. 
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L s k Ng connected disconnected glucball 

configs configs configs 



12 


0.1390 


24 


151 


301 


301 


1(5 


0.1390 


200 


90 


94 


390 


12 


0.1395 


24 


121 


253 


505 


16 


0.1395 


48 


100 


106 


424 


12 


0.1398 


24 


98 


169 


170 


16 


0.1398 


48 


69 


75 


298 


12 


0.13843 


24 


482 


100 


100 


12 


0.14077 


24 


482 


100 


100 



Table 2: Statistics of connected, disconnected and glueball calculations. For 
the fermionic disconnected correlation, the variance reduction methods with 
Ns samples were different for L s = 12 and 16 as described in the appendix. 



4 Scalar mesons 

Within the quenched approximation, there will be two distinct types of scalar 
meson: qq mesons and scalar glueballs. In full QCD, these two types of state 
will mix, resulting in the observed experimental spectrum of scalar mesons. As 
an aid to disentangling this experimental situation, we here explore the lattice 
predictions for scalar mesons. As a first step we evaluate the mixing matrix 
elements in the quenched approximation. This has been explored Q previously 
and here we discuss the problems associated with determining such hadronic 
mixing on the lattice. 

4.1 Quenched lattice results 

We explored this mixing in the quenched approximation (see Tables 1 and 2) 
using SW-clover valence quarks of two different masses Q . The zero-momentum 
glueball operators were measured at every time-slice in the usual way [[lof 
and the disconnected quark loops are measured as described in the appendix, 
namely with sufficient stochastic samples that no significant error arises from 
the stochastic algorithm. The connected quark correlators were taken from pre- 
vious measurements ^| . Since the scalar meson or glueball has vacuum quantum 
numbers, we subtract the vacuum contribution in the other types of correlation 
we measure. Our results for all of these types of correlation are illustrated in 
fig. | for the case of one choice of glueball operator and one (local) mesonic 
operator at our lighter quark mass. 

From the connected quark correlator C ss , at hopping parameters k — 0.14077 
and 0.13843, we find a scalar qq masses of 1.39(5), 1.36(2) respectively in lattice 
units, fitting local and two types of smeared operator to one state in the t 
range 2 to 8. Note that the mass ordering is not as expected (namely meson 
with lighter quarks being lighter) but the errors are large enough to cover near 
equality. This meson mass value is somewhat larger than that reported || at 
the same /3 value but using Wilson quarks of mass corresponding to strange (ie 
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Figure 2: Quenched scalar correlations with quark masses approximately 
strange. Here HH is the connected mesonic correlation (C ss ), GG is the glueball 
correlation (C mm ), DD is the disconnected mesonic correlation (D) and GD is 
the cross correlation between glueball and meson operators (H) . Lattice results 
are illustrated for one glueball operator and one (local) meson operator. The 
curves are a simple mixing model, as described in the text. 



our k — 0.14077), namely 1.29(2). This discrepancy is not surprising since the 
SW-clover formalism we use has improved control of order a effects compared 
to the Wilson discretisation. 

For the glueball mass, which is of course independent of fermion formalism in 
quenched studies, we use the higher statistics result |j| of 0.95(2) in lattice units 
obtained for t > 2. Our result for the glueball correlator (C mm (t)) is consistent 
with a single exponential with this mass for t > 1. Since our glueball correlator 
has large errors for t > 1, in evaluating the expressions shown in figs. || and 
[| we use our measured glueball correlations at t = 1 but for t > 1 we assume 
the glueball correlation has the mass dependence given by ma — 0.95 as found 
in higher statistics studies. As shown in fig. ||, this glueball mass lies below 
the continuum extrapolation because of order a 2 lattice artefacts. To convert 



12 



0.4 



0.2 



K=0. 14077 GD X 

DD + 



0.0 



1 2 
t 



Figure 3: The mixing coefficient xqo is determined in lattice units with a rs 1 
GeV, at each t value, from quenched scalar correlations with quark masses 
approximately strange: DD is the disconnected mesonic correlation (D) and 
GD is the cross correlation between glueball and meson operators (H). 



to physical units, we use r^ja = 2.94 and then conventionally tq « 0.5 fm, so 
a- 1 a 1.1 GeV. 

Then given these mass values, one can attempt to describe the disconnected 
and cross correlations (D and H) in terms of the one free parameter, the mix- 
ing strength xqq. The errors on our determinations of these correlations are 
quite large: 25% for D and 50% for H at t = 3. We use local fermionic op- 
erators in these comparisons since the smearings used in the determination of 
the connected and disconnected fermionic correlator were different for historical 
reasons. 

The curves shown in fig. ^ are from the lattice model described above with 
ground state contribution only and with xooa = 0.3. This is seen to give a 
reasonable overall description for 1 > t > 3. Alternatively, in fig. ||, || we give 
the value of the mixing parameter xqq for each t- value obtained from taking the 
data on H and D respectively and assuming the lattice mixing expressions of 
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0.4 



0.2 



0.0 



K=0. 13843 



GD X 
DD + 



Figure 4: As fig. || but with quenched scalar correlations with quark masses 
approximately twice strange. 



the previous section (eqs 32 



3|) with no excited state contributions. We find 
for each quark mass that £00 a ~ 0.3. 

Since our results for the diagonal correlations of glueball and meson operators 
(C rnm and C ss ) are only reasonably approximated by the ground state for t g > 1 
and t q > 1, as discussed above, the mixing should be studied from correlations 
with t 2. As shown in fig. we find that the signal becomes very noisy by 
t = 3 already. With 100 times the statistics we would be able to determine the 
t = 4 mixing correlations (D and H) to 10%. This implies that much larger 
data sets (number of gauge configurations about 10000) are needed to give a 
more definitive answer to the mixing in this case. Even then since there is only 
a power suppression of excited states, one would need precise data to large t to 
have small systematic errors from excited state contributions. 

The consistency of the determinations of x$o from different t values and 
different quantities does, however, act as a cross check that our results are 
consistent with the assumption that a single ground state dominates. Because 
of the lack of control of excited state contributions, we can only quote the 
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statistical error on the mixing xoo- Assuming no excited state contamination 
and taking t > 2, we obtain xooa = 0.26(4) for strange quarks, corresponding 
to xooro — 0.76(12). This is the quenched result for one flavour and we see no 
sign of any significant difference in xqo as we vary the quark mass since we have 
Xooa = 0.32(4) for heavier quarks. 

A previous work Q has studied glueball mixing with a scalar meson in the 
quenched approximation. They used several (3 values, Wilson fermions and 
concentrated on the cross correlation H to determine Xoo- At (3 = 5.7 and for 
quark mass near strange, we can make a direct comparison, bearing in mind that 
the order a corrections are significant at such a coarse lattice spacing and will be 
different for Wilson and SW-clover fermion formalisms, indeed the SW-clover 
formalism we use is focussed on removing these order a effects. 

At (3 = 5.7 they have t q — 2 and t g = 2 and they determine their mixing 
coefficient from data for H with t > 2 by assuming no excited state contributions 
to H. They do not consider data on D. Note that, as discussed above, from 
measurement of H alone, it is impossible in principle to confirm that excited 
state contributions are absent. Their quoted result for strange quarks is xooa = 
0.211(16). This is broadly compatible with our estimate of xooa — 0.26(4) 
bearing in mind that different fermion discretisations were used. 

At larger (3 (up to 6.2) their results for the mixing are that, at the strange 
quark mass, the mixing tends to a very small value in the continuum limit. The 
situation concerning excited state contamination is even worse at larger j3 since 
they find t q = 6 and t g = 4 at (3 = 6.2, whereas they still fit H for t > 2. Thus 
their mixing estimates at larger (3 are even more susceptible to excited state 
contamination. Moreover they do not make use of the disconnected correlator 
D to constrain their assumptions further. 

We find a significant mixing at coarse lattice spacings using a fermion formal- 
ism that has been shown to have reduced order a corrections. This implies that 
we would expect a substantial mixing in the continuum if order a 2 corrections 
were also to be small. This is in contrast to the conclusion [0 that the a depen- 
dence of the mixing (in physical units) is such that the continuum mixing is very 
small. We conclude that we find evidence for a mixing strength xoo^o — 0.76(12) 
with one flavour of quarks of strange mass at our lattice spacing. 

4.2 Full QCD scalar mesons 

To compare with our results using dynamical sea quarks with Nf = 2, we 
estimate the effects of our quenched determination of the mixing if applied to 
that case. For this estimate we take xooa = 0.3 for one quark flavour. Then for 
quarks of strange mass our unmixed states with moa = 0.95 and soa — 1.39 will 
be mixed by off diagonal element xooa = 0.3y/Nf =0.42 giving a mass mixing 
it 1 <\ t r * 

/ amo V2axoo \ ( 0.95 0.42 \ 
^ V2ax 00 aso )~\ 0.42 1.39 J 
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which gives mass eigenstates pushed apart to ma — 0.69 and 1.65 (and for the 
heavier quarks with the same mixing strength then 0.95 and 1.36 will be mixed 
to 0.68 and 1.63). Thus for strange quarks, the lightest scalar meson would be 
reduced in mass by approximately 0.24 in lattice units. If our quenched mixing 
strength were be to applied to the scalar mass matrix, it results in a downward 
shift for Nf = 2 of the lattice glueball mass by 25%. 

We now consider the sea-quark case explicitly, where the mixing will be 
observed directly from the resulting mass values. What we can measure in that 
case is the non-singlet mass and the ground and excited states in the flavour 
singlet sector. Based on the results from the quenched study, we would expect 
in the flavour singlet sector that the ground state mass lies considerable below 
the flavour non-singlet mass and the first excited state mass is slightly above 
the flavour non-singlet mass. 

We use 4 scalar meson operators (i) closed Wilson loops (glueball operators) 
of two different sizes (Teper-smeared 1 10 ) and (ii) qq operators which are local 
and separated by fuzzed links. For the fermionic correlations, we include the 
connected and disconnected contributions as given by the Wick formalism. We 
evaluate the 4x4 matrix of correlations. Since different numbers of gauge 
configurations are analysed for different operators, where necessary, we average 
results from nearby gauge configurations to have a consistent bootstrap sample 
for error analysis. 

We should also consider two particle states (two pseudoscalar mesons with 
angular momentum 1, for example) which can mix with scalar mesons and 
glueballs. On a lattice the lightest such state with overall momentum zero will 
have energy 2a(m^. + 4tt 2 / L 2 ) 1 / 2 which is at 1.53 for the case of L s = 12 and 
k = 0.1395 - this will apply to the flavour singlet case. For the flavour non- 
singlet scalar meson then the nr/ mode will be the lightest and this will be even 
heavier. These two particle energies are sufficiently heavy that we shall ignore 
these states in our present work. They will, however, become important as the 
quark mass is reduced and the lattice size is increased. 

The glueball and fermionic singlet correlations have a vacuum contribution. 
For the glueball fits we deal with this by using 2 state fits with one state con- 
strained to have mass=0 (namely the vacuum). In such fits it is important to 
use correlated fits to have a meaningful expression for \ 2 - Since we are fitting 
to many (up to 60) different types of data (different operators at source/sink 
and different t values), we use the technique of retaining exactly the N e largest 
eigenvalues of the correlation among the data set and setting the remaining 
eigenvalues equal |Tl| . This avoids spurious correlations being induced because 
of our limited sample size. We use N e = 10, 8 for L s = 12, 16 respectively. 
For our fits to the full 4x4 matrix of singlet correlations, since the number 
of gauge samples is different for different observables which makes the vacuum 
contribution depend on the observable, it is preferable to subtract the vacuum 
contribution and fit the resulting connected correlation. The errors quoted on 
the fits are statistical from bootstrap analysis and do not include systematic 
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L a 


K 


m GB a 


mpga 


m NS a 


12 


0.1390 


0.40(6) 


0.54(3) 


1.23(4) 


1(5 


0.1390 


0.53(7) 


0.47(3) 


1.19(5) 


12 


0.1395 


0.49(4) 


0.46(2) 


1.23(4) 


16 


0.1395 


0.70(9) 


0.75(4) 


1.18(8) 


12 


0.1398 


0.48(10) 


0.47(3) 


1.00(5) 


16 


0.1398 


0.58(8) 


0.66(4) 


0.99(6) 



Table 3: Ground state scalar masses from fits to the glueball sector (GB); the 
whole flavour non-singlet sector (FS: glueball and fermionic operators) and to 
the fermionic non-singlet sector (NS). The errors quoted are statistical only. 

errors from varying the fit range or fit function: these are at least comparable 
in size. 

We can now measure directly the scalar spectrum for Nf = 2 to explore 
this. Results are shown in Table [| We obtain the favour non-singlet mass 
(wijvs) from a two state fit from t = 2 to 7 to the 2x2 matrix of connected 
fermionic correlators. For the singlet, we now use both glueball and qq operators 
(vacuum subracted) and find an acceptable fit with 1 state to the 4x4 matrix 
of correlations for t from 2 to 7: the results are shown as mps® m 

Table § 

We find that the mass obtained from fitting only the 2x2 matrix of glueball 
correlations (mcs) is consistent with the full fit, as it should be. Moreover, we 
see a surprisingly low scalar mass - as emphasised in fig. ^| which compares with 
quenched results and the SESAM Nf = 2 values [pf. 

It would be interesting, as discussed above, to obtain the excited state mass 
in the flavour singlet sector. We expect this above toats at around am' = 1.4. 
We have used the variational method for t = 1,2 to extract the two lightest 
mass eigenstates from our 4x4 matrix of vacuum-subracted correlators. This 
variational ground state mass (0.44(1) for k = 0.1395) agrees quite well with the 
fitted value shown in Table ||, as expected. The next state is poorly determined 
although for the case with best statistics (L s = 12 and k = 0.1395) we find 
am' = 1.27 which is close to am^s- It will be interesting to explore this further 
with higher statistics. Note that it is difficult to determine this mass since the 
signal is swamped by that of two lighter states (the vacuum at m = and the 
ground state at am ss 0.5). 

We do expect a relatively light flavour-singlet scalar mass because of mixing 
effects as described above which would reduce the mass by 25%. This could 
explain in part our low scalar mass but other explanations are also worth ex- 
ploring. For example the order a 2 corrections might be anomalously large for 
our lattice implementation (e.g. twice as large as in the quenched Wilson case). 

Another possible explanation of the light flavour-singlet scalar mass we find 
would be a partial restoration of chiral symmetry in a finite volume, however, 
this would be expected only for very light quark masses. We do find that our 
flavour-singlet scalar meson masses are lighter than the pseudoscalar non-singlet 
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(pion) mass for L s = 12 for the range of sea quark masses considered here. The 
spatial size with L s = 12 is 1.7 fm and no evidence of finite size effects was seen 
here in a study of flavour non-singlet correlators . Indeed we see no significant 
sign of any spatial size dependence in the non-singlet scalar masses reported in 
Table |. 

To explore this further, we have made a study of flavour singlet correlators 
on 16 3 spatial lattices to check directly for finite size effects and the results are 
presented in Table ||. The signal to noise from zero momentum correlations is 
worse for the larger volume for glucball and disconnected correlations. Also we 
find that the excited state contributions are relatively stronger for our operators. 
Thus the systematic fit errors for L s — 16 are also considerably larger than those 
for the smaller volume. Even though the signal from the larger spatial volume 
is relatively poorly determined, we do see some evidence (at the 2a level if 
the sytematic errors are taken as comparable to the statistical ones) of a higher 
scalar mass (for mgg and tufs) on the larger spatial lattice. In order to explore 
larger L s values, the signal to noise can be improved by considering non-zero 
momentum correlators, as is the case for glueball studies pi. 

One conclusion is that it would be valuable to use a finer lattice spacing or 
an improved gauge discretisation so that any suppression of the glueball mass 
by order a 2 effects would be reduced. This would increase the glueball mass 
and hence reduce the magnitude of the signal we sec, but it would move the 
parameters into a region closer to experiment. 



5 Conclusion 

Hadronic mixing as exemplified by the glueball mixing with a scalar meson can 
be explored using lattice methods. In the quenched approximation, one can 
determine the mixing strengths although the systematic errors in this determi- 
nation are large as we have discussed. In studies with sea quarks, the mixed 
spectrum itself is obtained which gives complementary information. 

In a preliminary study of this glueball mixing, we find a large mixing in a 
quenched study and, consistently, a large suppression of the mass of the lightest 
scalar meson when sea quarks are included. These studies are computationally 
difficult and we have used a coarse lattice spacing, albeit with an improved 
(SW-clover) fermion formalism. It will be necessary to extend these studies to 
smaller lattice spacing in order to have more confidence in their relevance to the 
phcnomcnological situation. 
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0.00 0.05 0.10 

(a/r ) 2 

Figure 5: The scalar mass versus a . The quenched results p], fuj [D| @ 
are for the scalar glueball and are shown by boxes. The results from Nf = 2 
flavours of sea quark are from glueballs |l5| (crosses from SESAM) and the 
lightest flavour singlet scalar we find here (circles). 

6 Appendix: Disconnected fermion loop evalu- 
ation 

To evaluate disconnected quark loops with zero momentum, we need to sum 
over propagators from sources at each spatial location at a given time slice. This 
problem has been approached using stochastic source methods |l^, [l9| ^Cj . 
Here we describe in more detail the variance reduction techniques described 
previously For a similar method see ref. p2[ . 

To study flavour singlet mesons, we need to consider quark loops which are 
disconnected (often called hairpins), namely evaluate TrrM -1 where the sum 
(trace) is over all space (for zero momentum) at a given time value and all 
colours and spins. Here M is the lattice fermion matrix, T is a combination 
of the appropriate 7-matrix and a product of spatial gauge links if a non-local 
(fuzzed) operator is used for the meson. 
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Using a random volume source £ (where (£*£) = 1 for the same colour, 
spin and space-time component and zero otherwise) then solving Mcf> — £, one 
can evaluate unbiased estimates of the propagator M~y from (£y4> x ) where the 
average is over the Ns samples of the stochastic source. The computational 
overhead of this method lies entirely in the inversion of M to obtain <j> from £ 
for each of the N$ samples. 

The drawback of this approach is that the variance on these estimates can 
be very large, so that typically hundreds of samples are needed. Here we present 
a method which succeeds in reducing this variance substantially at rather small 
computational expense. 

The variance reduction is based on expressing the fermion matrix M as 

M = C - D = C(\-C~ 1 D) = (l-DC-^C (34) 

where C is easy to invert, for example the SW-clover term which is local in 
space or in the Wilson case where one can chose (7=1. Then we have the exact 
identity 

M~ l = C- 1 + C^DC'- 1 + ... + (C^D^C- 1 

+ (C- 1 L>)™ 1 Af- 1 (L>C- 1 )™ 2 (35) 

with n± + ri2 = n = m + 1. 

Our strategy will be to evaluate the last term in this expression stochastically 
and to evaluate the preceding terms as exactly as possible. We will refer to these 
terms in our following discussion as the stochastic and exact terms. If these 
exact terms can be evaluated precisely, then it is plausible that the stochastic 
term will contribute less variance to the overall estimate of M _1 than in the 
m = case where there would be no such exact terms. These exact terms can 
be evaluated either directly (for example terms with odd powers of D vanish 
in the evaluation of a local trace) or as a subsidiary stochastic calculation with 
more samples to achieve good precision and at relatively small computational 
overhead since no inversion is required. 

Since this approach is a variant of the hopping parameter expansion, it 
might be suspected that the convergence was poor since at each higher order 
8 extra terms are present with coefficients which are of order k w 1/8. In our 
application, however, these 8 terms contribute with random strengths - like a 
random walk. So they have an effective weight which is more like y/E which is 
smaller. So we do find a reduced variance including more terms, but at the cost 
of some extra computation in evaluating the additional exact terms on the right 
hand side of eq. |B[ 

A special case of this (ni = nz = 2) with Wilson fermions (for which (7=1 
and the terms with up to 3 powers of D vanish for TrM -1 ) employing Gaussian 
noise was used by the bermion group Jl7[ previously. 

Using the stochastic volume source, the variance reduced expression can be 
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rewritten (assuming C is hermitian) as 



r-l 

yx 



... + {{(C- 1 D^OlT xy ((C-'D) n ^) v ) 



(36) 



X 



so the stochastic term may be evaluated as an average over stochastic samples £ 
after inversion to obtain cf> = M^ 1 ^. In the application of this paper we take T 
to have the Dirac structure / whereas in other applications jl6], we consider 
75 and 7574 also. 

Only the even exact terms in the series when T is local are non-zero and we 
calculate the m = and 2 cases explicitly. For hadron operators with fuzzed 
paths of length np the series starts at m — np (this we calculate explicitly) 
and then has alternate terms zero. The explicit calculations referred to are 
rather cumbersome for clover fermions, so in some cases we actually evaluate 
the simpler Wilson expression and then evaluate the difference stochastically. 
The generic non-zero terms (A) in the series were calculated stochastically using 
ANs Z2 noise samples £ using 



Taking n\ ps n% gives an averaging over a smaller volume than taking an 
asymmetric choice. We find that an asymmetric choice gives a smaller variance, 
presumably because it does involve averaging over a larger volume. For different 
disconnected observables, the optimum strategy is not necessarily the same. In 
this work, we find a overall good choice to be n\ = 0, n-i = 16 with Gaussian 
noise. This results in a reduction by a factor of 0.2 to 0.3 of the standard 
deviation of the samples. 

Using larger values of n\ and ni implies that the estimate of M~ x is very 
non-local, involving £ values up to n\ and rii away. To evaluate correlators 
between traces at t\ and £2, one must require that the samples of stochastic 
volume source used in the two cases are different so that there is no bias. We 
use N$ — 24 stochastic samples and this condition is readily implemented with 
essentially no loss of statistics. This number of samples was chosen to make the 
stochastic sampling error smaller than the intrinsic variance from one time slice 
(for example at L s — 12 with k = 0.1395, the ratio of the stochastic sampling 
error to the standard deviation over time slices was 0.5, 1.0, 0.2 for T = 75, 
7574, I respectively for local operators and the ratio was about 50% bigger for 
fuzzed hadronic operators). The computational effort in terms of the number of 
inversions is equivalent to that in obtaining two conventional propagators (from 
two sources of all colour-spins) . We also make use of the fact [|o) that one can 
truncate the inversion at a larger residual without measurable bias, using 10~ 9 
in an MR inverter as residual/source. 

In conclusion, using N$ — 24 samples for each gauge configuration corre- 
sponds to 2 conventional propagator determinations (from all 12 colour-spin 




(37) 
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combinations) and so is not a particularly big computational challenge and yet 
the resulting measurement of the disconnected fermion loop has a stochastic 
error which is unbiased and less than its intrinsic error. 

For spatial size 16, we made use of existing solver codes and chose not to 
use the variance reduction method described above, using instead m = n 2 = 0. 
To achieve some variance reduction, we used the method |24|] of using pairs 
of sources with the same random Z2 numbers but with the second of the pair 
multiplied by 7,75 where i is chosen randomly from 1 ... 3. This has the effect 
of reducing the standard error by a factor of two for r = 75 for only a doubling 
of CPU effort. In this case we used either N$ — 48 or 200, the larger value 
being motivated by the need to get a more precise estimate for T — 7574. 
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